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Abstract 

We reformulate the standard local equations of general relativity for 
asymptotically flat spacetimes in terms of two non-local quantities, the holon¬ 
omy H around certain closed null loops on characteristic surfaces and the light 
cone cut function Z, which describes the intersection of the future null cones 
from arbitrary spacetime points, with future null infinity. 

We obtain a set of differential equations for H and Z equivalent to the 
vacuum Einstein equations. By finding an algebraic relation between H and 
Z this set of equations is reduced to just two coupled equations: an integro- 
differential equation for Z which yields the conformal structure of the under¬ 
lying spacetime and a linear differential equation for the “vacuum” confor¬ 
mal factor. These equations, which apply to all vacuum asymptotically flat 
spacetimes, are however lengthy and complicated and we do not yet know 
of any solution generating technique. They nevertheless are amenable to an 
attractive perturbative scheme which has Minkowski space as a zeroth order 
solution. 
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I. INTRODUCTION. 


The vacuum Maxwell fields on Minkowski space can be expressed as Kirchhoff integrals 
taken over initial data, i.e., data given on a Cauchy surface. If we replace the spacclikc 
surface by a characteristic, or null, surface, we get the analog to the Kirchhoff integral, the 
D’Adhemar integral [Q. Maxwell and Yang-Mills fields on both flat and asymptotically flat 
spacetimes have been given such D’Adhemar-like reformulations with considerable success 
The present work is an attempt to extend this technique to general relativity. We 
derive differential equations whose solutions would allow us to express the gravitational field 
at any given spacetime point as a series of D’Adhemar-like integral over free data given on 
a characteristic surface. 

First consider the well understood analog to our problem, determining the Maxwell field 
at some interior point on Minkowski space. Free data is given on the null surface X, which is 
the three dimensional boundary of the spacetime manifold Ad. The future light cone ema¬ 
nating from a given interior point x a would intersect the future half Z + of X, the intersection 
being a 2-surface. We refer to this 2-surface as the “light cone cut” of X + associated with 
x a . The D’Adhemar formulation of the Maxwell equations involves integrals over the free 
data given on the light cone cut. It is crucial that the light cone structure of the underly¬ 
ing manifold be known, in order to perform this integration. Although for the Minkowski 
metric this is readily available, it takes considerable effort to obtain this information for an 
arbitrary background. A generalization of the D’Adhemar formulation of the Maxwell field 
in Minkowski space to the case when the background is curved is available ||, although the 
resulting equations become complicated. A further generalization of this technique to Yang- 
Mills fields on both flat and asymptotically flat spacetimes are also available ||. In these 
(Yang-Mills) cases, in addition to the above mentioned difficulty in the determination of the 
light cones, the final equations are further complicated by the introduction of non-linearity. 
The D’Adhemar integrals (for the Yang-Mills case and the Maxwell case on a curved back¬ 
ground) shows explicitly that propagation becomes non-Huygens-i.e., propagation is not 


2 


confined to the light cone. 

A D’Adhemar-like formulation for general relativity might seem inaccessible since the 
gravitational field is not separable from the background and its characteristic surfaces: it 
requires simultaneously solving for the field and the null surfaces. It is nevertheless possible 
to give such a formulation by dividing the equations into two parts: one resembling the 
D’Adhemar integrals for a (self-dual or anti-self-dual) Yang-Mills field, and the other yielding 
the null surfaces that are needed to give meaning to the integrals in the first part. The two 
parts are strongly coupled. 

The first set of equations is equivalent to the D’Adhemar version of the self-dual and the 
anti-self- dual Yang-Mills equations for the self-dual and anti-self-dual parts of an 0(3,1) 
connection, respectively, on a spacetime manifold A4 with an arbitrary asymptotically flat 
metric. The basic idea is to treat a special case of the 0(3,1) Yang-Mills connection so that 
it agrees with that of the background gravitational connection. This is done by imposing 
very restrictive conditions on the Yang-Mills data. 

The second set of equations determines the light cones of the manifold A4. For a given 
metric, one could have obtained this information by integrating the null geodesic equations 
(see for example |7J]). However, in our case where the metric is as yet undetermined, we 
derive differential equations for the characteristic surfaces that are coupled to the 0(3,1) 
Yang-Mills field equations. 

The above two parts, which look hopelessly intertwined at first, nevertheless give results 
that are tractable and comprehensible. The exact equations, while complicated, are in 
principle solvable. An attractive feature of this approach is that the equations are amenable 
to an approximation scheme wherein the n-th order correction to the field is expressible as 
D’Adhemar integrals in terms of lower order. 

As the material of this paper is not at all close to mainstream ideas, we believe that it 
might be appropriate first to give a broad perspective before going into the details. 

There are two different sets of ideas that we try to weld together. The first is based on 
the use of characteristic surfaces on arbitrary asymptotically flat Lorentzian spacetimes. We 
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make essential use of the future light cones N x from arbitrary spacetime points x a and their 
intersections, C x , with future null infinity, Z + , (referred to as the light cone cuts)-as well as 
the past light cones from points on Z + . These characteristic surfaces are used to give the 
D’Adhemar-like reformulation of the Yang-Mills equations. Using Bondi coordinates (u, (, () 
on 1 + ,C X is described by a function referred to as the light cone cut function on X + of the 
form u = Z(x a parametrized by the interior points x a (see Figures 1 and 2). The cut 
function plays a dual role: for a fixed interior point x a it describes the “cut” of Z + , and for a 
fixed point ( u , () on J + , it describes the past light cone of that point as x a is varied. From 

the latter point of view the cut function becomes a three-parameter family of characteristic 
surfaces. The conformal structure of spacetime can be recovered from knowledge of the cut 
function, i.e., the conformal metric can be explicitly written in terms of derivatives of Z. 
More specifically, the conformal metric can be completely expressed in terms of A, defined 
from Z by A = d 2 Z. The quantity A, becomes one of our basic variables. These ideas are 
developed in Section II. 

The second set of ideas of concern to us (see Section III) is based on the fact that, in some 
sense, the 0(3,1) Yang-Mills equations contain, as a special class of solutions, the vacuum 
Einstein equations. The full Yang-Mills curvature, F, can be broken into four pieces; self¬ 
dual and anti-self-dual on the internal indices as well as self-dual and anti-self-dual on the 
spacetime indices, so that 

F = + F + + + F~ + ~F + + ~F~, ( 1 ) 


where the +, — symbols to the left refer to spacetime duality and those to the right refer to 
internal duality. The Einstein equations are almost encoded (Tip into the algebraic statement 
that + F~ = ~F + = 0 needing only an appropriate restriction of the data. The basic idea 
is to think of solving these equations (for the self-dual and the anti-self-dual fields + F + and 
~ F~ and their respective connections) on a given asymptotically flat vacuum background 
metric, then introduce a soldering form and restrict the data so that it coincides with the 
background connection data. An important technical device is to give these equations a 
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D’Adhemar-likc formulation which, as was mentioned earlier, can be done for Yang-Mills 
fields. The tool for this is the use of the holonomy operator H (around special null paths) as 
the primitive variable for the Yang-Mills fields. The Bianchi identities, equivalent (because 
of the self- and anti-self-duality) to the Yang-Mills held equations become one set of our 
final equations. They are symbolically written as equations for H as 

B g (H,A) = 0, (2) 

with A being the restricted characteristic data. 

When the soldering form (a tetrad held) is introduced, we can obtain a set of relationships 
between the spacetime variable A = d 2 Z, and H , that we refer to as the A (H )-relations given 
symbolically by 

A = A (H). (3) 

Our third set of equations, referred to as the “held equations”, are simply the algebraic 
equations 

+F~ = ~F + = 0, (4) 

expressed as functions of the components of the holonomy operator, symbolically written as 

V(H) = 0. (5) 

Our final task is two-fold: To 

1. eliminate the holonomy variables from these three relationships and be left with equa¬ 
tions only for the characteristic surfaces, i.e., equations for Z (Section lll.E.) These 
equations constitute the conformal Einstein equations and yield an equation that we 
refer to as the light cone cut equation (LCCE) which involves only A and the free data, 
and 

2. find an equation, derivable from equation (|j), for the conformal factor that converts 
the conformal metric to a vacuum metric. 

In Section IV we discuss the linearization of the theory and a perturbation scheme. 
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II. LIGHT CONE CUTS. 


Let us consider a real, four-dimensional, asymptotically flat spacetime, M, with a con¬ 
formal metric g ab and future null boundary X + . We use Bondi coordinates (w, C,C) to 
coordinatize X + . A coordinate grid of the (u =const.) and the (£, ( = const.) surfaces on 
X + is shown in Figure 1. 

Consider the future light cone N x emanating from an internal point x a . A relation among 
the coordinates such as u = /((,() describes locally a 2-surface which we refer to 

as a ‘cut’. The intersection of N x with X + is a special cut, the ‘light cone cut’ C x (Figure 
2), given by 

u = Z(x a , C,0- (6) 

We refer to Z(x a ,(X) as the light cone cut function. The function Z(x a ,(,() is our 
basic variable. The conformal metric can be given as an explicit function of Z 0. 

One could also interpret equation (;6|) as a description of the past light cone of the point 
(w,C,0 cm Z + (also shown in Figure 2). That is, keeping u, ( and ( fixed, if we vary x a , 
this equation is the locus of all spacetime points x a that are null connected to the point 
(w, C,C), which by definition is the past light cone of °n X + . From this important 

observation, we see that for fixed values of (u, C,C)> V a Z is a null covector, i.e., 

g ab z, a (x,c,C)z, b (x,c,0 = o. (r) 

Thus it follows that Z^ a (x, (,() at the point x a sweeps out the null cone at that point 
as we vary ( and (. Though it has been discussed elsewhere 0 and we will not go into the 
details, it is from equation 0), by taking several ( and ( derivatives, the entire conformal 
metric can be reconstructed completely in terms of Z or more specifically A = 5 2 Z. 

The sphere of null directions at the point x a is coordinatized by ( and ( . We use these 
coordinates for S 2 instead of the usual (9, (ft) because covariant differentiation on the sphere, 
which appears in many of the equations, takes on a particularly simple form in terms of 
(C,0- We introduce the operators 3 and 8 0 
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( 8 ) 


Sa.Hl + CC) 1 S -^~ [(i + CC) S « S ] 


and 


8a s = (1 + CC) 


l+s 


d_ r 

dC 


(i + CC) a s 


(9) 


which operate on the quantity a s defined on the sphere, where s, the spin weight, is assigned 
according to how a transforms under a specific transformation. (See 13 for more about 
this transformation and properties of 8 and 8.) It helps to think of 3 and 8 loosely as 
differentiation with respect to C and C respectively. 

Assuming that we know the light cone cut function, and with the above definitions of 3 
and 8, the following set of quantities are well-defined: 


u = Z(x a , C,C), 

W = 5Z(x Q ,C,C), 
uJ=BZ(x a X,(), 
R = mZ(x a XX), 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


with (u,R,lv,u) having, respectively, spin-weights (0,0,1,—1). This set defines a sphere's 
worth of coordinate transformations on the spacetime parametrized by (C,C). Let 

(9°,9 1 1 9 + 1 9~) = (u,R,u,uf). (14) 

With this notation, the (C, C)-dependent coordinate transformation can be written as 

9 l — 9 i (x a , C, C)i (15) 

where the indices i,j will take on the values {0,1, +, — 

We now construct 

d 2 z = \(x a xX). (16) 


The quantity A and its conjugate can be expressed as functions of the 9 l by inverting 
the transformation (jl5]) and eliminating the x a , to obtain 
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A(^C,0 = 9 2 ^. 


Later we show that the (conformal) Einstein equations can be encoded into the choice 


With the 9 l coordinate system, we have a one-form basis: 


(W l = 9\dx a , with 9\ a = 9\ 


and the dual basis 


d9 i 1 dx a ’ 


with the relations 


9\9 a j = <W, 


9\9 b i = S a b . 


Any vector V a can be written as 14 = ld$4 and in particular we have 


A, a — Aj 04 — A fiZ a + etc... 


The A ;i play a basic role in what follows. 

Note also that 3 9 l a can also be expressed as a linear combination of the 9 l a , i.e., as 


3 9\ = T l j9 3 a 


From 9\ = {Z a , ddZ a , dZ a , dZ a }, we obtain, e.g., that 


3 9° a = 9 + a , 

3 9 + a = A, 0 = A/4, 
3 9~ a = 9\, 


= T\9\ 


with 



lr 


T ■ = - 

- 1 - l 


1 - - A 1 (9A i + Ao<5 + , + K.5\ + A + Aj - 2<$“*) 

q L 

+8A, + A 0 <5“ i + A+5 1 * + A_ A ? ; - 2<S + i 


and 


q = (1 - AiAi). 


( 28 ) 


Likewise, 


se-„ = t 


(29) 


with T j obtained from the complex conjugate coefficients of T l j. The explicit form of these 
matrices with row and column indices {0,1, +, —} are: 

/ n n i n \ 


rpi _ 

1 3 ~ 


0 0 10 

T 1 o T\ T l + T 1 . 
A 0 Ai A + A_ 
0 10 0 


(30) 


and 


T ■ = 

1 3 


^ 0 0 0 1 ^ 
t\ t\ t\ t\ 

0 10 0 

y Aq Ai A + A_ J 


(31) 


III. HOLONOMY AND EINSTEIN EQUATIONS. 

A. 0(3, 1) Yang-Mills Equations. 

We begin with a real four-dimensional Lorentzian manifold JA where we assume that 
the light cone structure of the previous Section is known and study an 0(3,1) Yang-Mills 
held. The idea is to (roughly) think of the associated vector bundle as being the tangent 
bundle of the spacetime, though of course a soldering form is needed to make this precise. 
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The connection one-form, which is antisymmetric in the Lorentzian indices, i,j can 


be decomposed into its self-dual and anti-self-dual parts, 


7« J i = 7fi + 7: 




(32) 


where self-dual and anti-self-dual are defined by 


T„ ±f j = \ (7 Jj =F h’Jj) , 


(33) 


and duality by 


.*27 1 

7a = ^ 6 7a , 


(34) 


where is the alternating symbol with 60123 = — 1. The curvature tensor can be similarly 
decomposed as, 


ab j 


= F£j + F r 


ab 3"> 


(35) 


where the self-dual and anti-self-dual parts of F are constructed from the self-dual and 
anti-self-dual connections respectively, i.e., 

F abj = v \alf 3 + [7f, 7 6 ± ] V (36) 


With the above decompositions, the Bianchi identities and the field equations become 


and 



(37) 


+ ( 38 ) 

respectively. One is thus dealing with two independent Yang-Mills connections and fields. 

It is possible to further decompose each of the two curvature tensors, now on the space- 
time indices, into its spacetime self-dual and anti-self-dual parts, where we have used the 
existence of the Lorentzian metric. We will refer to spacetime duals as left duals and internal 
duals as right duals. The full curvature then has four terms: 
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1. the left and right self-dual part, + F+ b \ 

2. the left anti-self-dual and right self-dual part, ~F+ b \ 

3. the left self-dual and right anti-self-dual part, + F~ b , 

4. the left anti-self-dual and right anti-self-dual part, ~F~ b . 

Parts (1) and (2) are coupled as are parts (3) and (4), in the sense that they depend 
respectively on the 7 + and y - . 

We next assume that: 


~rr+ —+ rr~ — n 

r ab ~ r ab ~ U ’ 


(39) 


From this the field equations (|38|) are automatically satisfied via the Bianchi identities 
(ft7|). Thus we have two independent Yang-Mills fields: a (left) self-dual field, + F+ b and a 
(left) anti-self-dual field, ~F~ b satisfying the Bianchi identities (|37l) and the above condition 
(j39|). These two equations are rewritten in terms of a new variable, viz., the holonomy 
operator, in the next section. 

Equations (|37| ) and Q39] ) with a special choice of data (on Z + ) are equivalent to the 
Einstein equations. 


B. Holonomy and the Bianchi Identities. 

The variable that we have been (loosely) referring to as the holonomy operator and 
denoting by if, is more accurately the difference between the holonomy operator associated 
with an infinitesimal path and the identity. It should really be called the infinitesimal 
holonomy operator or the differential holonomy operator, though we will continue with the 
original name. There are two distinctly different sets of paths A x (£,£) and A x (£, £) (with 
their own holonomies FI and H) defined in the following manner: Consider an interior point 
x a of M.. For A x (£, £) we choose two null rays on the cone N x that are infinitesimally 
separated, namely, £ X (C,C) and Cz(C + extending from x a to Z + (see Figure 3). We 
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then form a closed loop by connecting the end points of these two rays on X + . At X + the 
two-form constructed from this connecting vector and the tangent vector to the geodesics is 
self-dual. (In Minkowski space the entire path lies in a self dual blade; it is the loss of this 
property in curved spacetimes that is the source of non-Huygens propagation for linear rest- 
mass zero fields). In a similar manner, one can choose the anti-self-dual triangle A x (£,£), 
which has £ x ((, C) and G(( + d(, () for its sides. 

For our 0(3,1) connection, the vectors that are being propagated around the closed loops 
(A x and A x ) are thought of (using a soldering form A a i(x) introduced below) as being in the 
tangent (or cotangent) bundle. The effect of applying the operators H or H to an arbitrary 
vector at X + is, respectively, 

V ,fl = V U (cFA + 



and 


V 


V v 



2 P 


(41) 


where the P = (1 + (() is for notational purposes. H takes us from a point on X + along 
£ x ((, C) down to x a then back to X + along + d(, Q and finally back to the starting point 
along a connecting vector on X + (see Figure 3). For parallel transport around this loop one 
obtains the following: 


H\(x‘XX) = A\+(G- 1 )*<iG\, 


(42) 


where G^ a (x a ,(X) is the parallel propagator that takes vectors from the point x a to X + 
along the null geodesic ^(C, C)> an d A'T are the asymptotic values of the components of the 
Christoffel symbol in the direction of the connecting vector on X + . See [|lTJ and [|l^] for a 
derivation of this result. 

It will be convenient (see (59)) to express the above relation in terms of a given null 
tetrad A a i(x) defined on M. and X + with normalization 


9ab A ?A j Viji 


(43) 
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with 


Vij = 


0 10 0 

10 0 0 

0 0 0 -1 

0 0-10 


\ 


(44) 


and H l j = H^ U X\X U j, etc. Tetrad indices are raised and lowered with rf ] and rjy, respec¬ 
tively. For the tetrad fields X^i at X + we will use the Bondi tetrad associated 

with the Bondi coordinates. 

The Bianchi Identities (37) can be written in terms of the holonomy operator as 


3 {H - Af s - 3 {H - Af 3 + [H — A, H — Af 3 - [H, A]* 3 = 0. (45) 


A derivation of this result, via an integration of (37) over the region V of Figure 3, can 
be found in (1J. Using to raise one of the indices, we obtain H'\ A l \ and their complex 
conjugates. Since these are skew in the ij indices, they can be separated into self-dual and 
anti-self-dual parts, yielding 

= and A iS = + AW (46) 


and 


H' s _ H Mti + h < )tj and A v = A (-)!} + A (+y,j 


(47) 


We have three complex non-trivial components of H^ 1 which we denote by H a and likewise 
three h a , where a = {1,2,3}, as shown in Table 1. As for the characteristic data, the A’s, 
most of them are zeroes and ones: they are given by the following: 


( 0 0 - 10 ^ 
0 0 1 — & B 

—& B 0 0 0 

v 1 -1 0 0 y 


( 48 ) 
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and 


/ 


*j = 


\ 


0 0 0 -1 

0 0 — W B 1 

1-10 0 
— a B 000 

The self-dual and anti-self-dual parts of H can be written out explicitly as well. 


V 


(49) 


) 


( 


tf(-)b = 


0 


\{H m -H+~) 


— ^(i4°i — H + ~ 
-H°+ 

0 


and 


h (m = 


-\(H m + H+- 
0 

-H°~ 


0 

0 

-H 1 - 

0 

-H l+ 

0 


H°+ 

0 

0 


1 (H 01 - H+- 

0 

H 1+ 

0 


0 

H 1 - 

0 

H 0+ 

0 

1 (H 01 + H+~) 
0 


\ 


(50) 


(51) 


and likewise for H. 

In order to raise and lower the i,j indices, we have to use the null-coordinate version of 
the Minkowski metric r/y and rfK For example, 


H °i = H 00 = 0, by skew symmetry , (52) 

H\ = -H°~. (53) 

Also note that the complex conjugate of a quantity, say, H°~ gives us H° + : the 0 and 
1 indices are insensitive to complex conjugation, while + and — are interchanged. The 
following table displays our notation. 
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new 

hJ 

spin 

new 

ij 

spin 

notation 

component 

weight 

notation 

component 

weight 

Hi 

-HH) o+ 

2 

Hi 

-#(+) o- 

-2 

h 2 

-HH) oi 

1 

h 2 

_jy(+) 01 

-1 

H 3 


0 

h 3 

-H(+n+ 

0 

hi 

_/ 7 ,(-)°+ 

0 

hi 


0 

h 2 


-1 

h 2 

_/ z (+) qi 

1 

h 3 

-h^n- 

-2 

h 

_/ 1 (+) i + 

2 


Table 1: Notation. 


With this notation the holonomy Bianchi identities fl45|) become: 

9/ii - QHi + 2h\H\ - 2H 1 h 2 + 2 H 2 = 0, (54) 

Bh 2 — BH 2 + h 3 Hi — H 3 (hi + 1) + hi = —u b Hi, (55) 

Bh 3 — BH 3 + 2h 2 ( K H 3 — 1) — 2H 2 h 3 = —Bu B + 2cr B H 2: (56) 

where we have the data or “driving” terms (i.e., those that involve a B ) on the right. Equation 
(|45|) actually contains this triplet as well as their conjugates. They are the result of taking 
the self-dual and anti-self-dual parts of equation (|45|). This is the first set of equations that 
was given symbolically as equation (||D in Section 1. (Note that the first equation is algebraic 
in h 2 and H 2 , and the second is algebraic in h 3 and H 3 . We will exploit this structure in 
Section III.D where we study these equations in more detail. 

C. The A(hf)-relations. 

In the previous Section we defined the parallel propagator, G^ a {x a , (), which takes 
vectors from a point on X + to an interior point x a along a null geodesic £ X (C, ()■ In particular, 
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the tetrad fields {A* M , A'fi} that are defined onI + can be parallclly propagated in this manner 
to the point x a . The parallclly propagated tetrad fields at x a are given by 


J — X* Qf* 

° a uwr a? 


(57) 


where A \ represents the tetrad {i^, n M , on Z + . We use the following notation for 

the parallclly propagated tetrad at x a : 


6 a — ai 6 a; 6 a; 6 a} — ^lai ^ai ■ 


(58) 


If we take two sets of tetrads on J + , one at (C,C) and the other (on the same cut, 
associated with x a ) at (C+dC, C)> and parallclly propagate both to the point x a , the difference 


in terms of the holonomy operator H and the free data A [11] is given by 


* — m ak 

a 11 3 j i 


(59) 


= H', - A‘ r (60) 

We have introduced three different tetrad bases: X\, X\, and e\. The first of these, X l a , 
is an arbitrary tetrad given at x a and therefore independent of ( and (. The second, A* M , is 
a ((, (^-dependent Bondi tetrad given on J + . And the third, e\, is the result of parallclly 
propagating the X\ in from X + to x a along £ X (C, (), and therefore is also ((, C)-dependent. 

From Section 111, we have at our disposal the Q l coordinates and the associated bases 9 l a 
and 6 a i at x a . Since any vector at x a can be expressed as a linear combination of either of the 
two sets {9 l a} and {e* a }, we can go from one to the other, using the invertible transformation 

e\ = EVV ( 61 ) 

We choose one of the ‘legs’ of the parallclly propagated tetrad vectors, e° a , equal to the 
vector V 0 Z, i.e., 

e° a = L = Z, a = 6\. (62) 
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Using this assumption and equations (^) and (|6(]) with our notation from Table 1, the 
transformation matrix = £ and its inverse £ l i = £ _1 can be calculated (by application 
of 5 and 8 to (59) and (60)) and written explicitly as functions of the hT s and h’s |T[. £ is 
given by 


/ 


PS) = 


1 0 0 
S x 0 S 1 ! SV 
■( H 2 + h 2 ) 0 (1 + hi) 


0 

Hi 




^ — {H 2 + ^ 2 ) 0 H\ (1 + hi) j 


(63) 


with 


S X o = -3 (H 2 + h 2 ) + (H 2 + h 2 )(H 2 + h 2 ) + H 1 (a B + h 3 ) + (1 + hi)(hh 3 - 1) 

S 1 ! = (1 + hi)(l + hi) + H 1 H 1 

EV = &Hi - (H-2 + h 2 )( 1 + hi) - Hi(H 2 - h 2 ) 

U 1 - = 5hi + (H 2 - h 2 )(l + hi) - Hi(H 2 + h 2 ) 


Using £, its conjugate, the equations (|23|), i.e., 

= Tj&a (64) 

and their conjugates, the T 1 / s and T*/s of Section II can be expressed directly in terms of 
the hTs and h’s via the by 

T = (8£) • S' 1 + £ ■ (H- A) • 5T 1 , (65) 


where the boldface represents the corresponding matrices [|IJ. Since T l j was shown in equa¬ 
tion ([31]) to have the form 


/ 


0 


0 


1 



T = (Tj) 


T x o T 1 ! T l + T 1 . 
Aq Ai A + A_ 


y 0 1 0 0 y 


( 66 ) 
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we obtain from (|65D and (66) the A(i/)-relations: 


Ai = Ai = 


2{h\ + 1)H\ 


{h\ + l)(hi + 1) + H\H\ 


A + = W - - (AiVh + 5Ai - 31ng) , 


A_ = - 
2 


3Ai 


A miH i) 
1 


Axd In 


1-V? 


A 0 = ^ + — ^S5Ai + 0{H 2 ), 


where W is dehned by 


T „ (3Ai — 2A_) 31 

W =-T-ding, 

Ai 


with 


<? = (! — AiAi). 


(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 


See [1] and |1 for a more complete derivation of these equations. The precise form of 
equation (711) turns out to be rather lengthy and is given in Appendix A of [|T|. 


D. The Holonomy Field Equations. 


In this section, we describe an important result obtained earlier |lT|, wherein the vacuum 
Einstein equations in the form 


+ f~ — n 

r ab ~ U 


(73) 


were imposed on the holonomy operator. The goal is to obtain held equations for the 
holonomy operator, and therefore we first need a relationship between the curvature tensor 
and the holonomy operator. 


It can been shown [|TT[ , using a non-Abelian version of Stokes’ theorem, that the holonomy 
operator satisfies 


H= / (F+ + *3 rM b is = fc(+> + H'■-> 


(74) 


'so 


18 


















and 


H = 


« + F-J rM b ds = i/<+> + fc<->, 


(75) 


'«0 


where the H and h with the pins and minus signs were defined in Section III.B, s is an 
affine parameter along the generators of the light cone at x a with s = s 0 at x a . M a and M a 


are connecting vectors between, respectively the long legs of A x and A^. See |11] for more 
details. 

The two equations ([74|) and ([7^) can be inverted to obtain F~ b and Ff~ b explicitly in terms 
of the components of the holonomy operator. Specifically F~ b and F+ b can be expressed as 
two derivatives of the H 1 s with respect to 9 1 = R. Using this expression for the curvature in 
+ F~ b = 0 one obtains a differential relationship between the different components of H and 
H, given by 


[q ha,R],R + S[q Ai h ai R\,R—[q A-iH a ,R],R + 8[q H^r]^, 


(76) 


where 

Ai 

Note that the equations ( 0 ) , which we refer to as the field, equations , are three linear 
relations between the H and h with coefficients depending Z via Ai. They are our final 
set of equations. In the next section, we will use these with the Bianchi identities and the 
A(/7)-relations to simplify the overall structure of the theory. 



E. The Full Theory. 

We have at this point three sets of equations: 

• the three holonomy Bianchi identities, 

Shi — QH\ + 2 h\H\ — 2Hih 2 + 2 H 2 = 0, 

Sh 2 - m 2 + hzH, - H 3 (h i + 1) + h = -F b H u 

Sh -3 — 3773 + 2h 2 (H 3 — 1) — 2H 2 h 3 = —Scr B + 2<r B H 2 , 
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(78) 

(79) 

(80) 






the holonomy field equations, 


[q h aj R] t R + 8[q Mha, R },R — [q A ,RH a:R ] tR + 5[q H a}R ] )R , 


(81) 


the A(i/)-relations„ 

A, = A, = — 


A_ = - 
2 


2 (hi + l)Hi 

(hi + l)(/ii + 1) + HiHi 

(82) 

^ (AiW + Mi-mnq) , 

(83) 

U-A 9 )j 

HiHi 

(84) 

jA, - jsSA, + 0(H' J ), 

(85) 


and their conjugates. The first set involves the H 's and /i’s and the data only, while the 
other two involve the H’s and K s and the A*. We can think of the FT s and h 's as describing 
the “field” and the A,, the “background”. In this Section we study the above equations and 
show how they can be reduced to a smaller and simpler set. 

This simplification procedure begins with the following important observation: equation 
(j32j) its complex conjugate are two algebraic relations between the six quantities, viz., Ai 
and Ai, Hi, Hi, hi and hi. Solving for hi and h\ algebraically yields, 

hi = - 5 ±- 1 , (86) 

o 

hi = ~~ 1, (87) 

where 


5 = '^-— K - and q = 1 — AjAi, 

Ai 


( 88 ) 


Thus hi and hi are not independent quantities and equations (|86|) and (^7|) allow us to 
eliminate them completely. This algebraic structure extends to the other h 1 s. We are able to 
solve, from equations (|78| ) and (j79j. for all the remaining h’s and completely eliminate them 
from the remainder of the analysis. The algebraic relationships between H a and h a for all 
a = {1, 2, 3} takes the form 
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( 89 ) 


u _ _ £U i r 

11 a. — i ^ oli 

0 

where G a is a function of Ai and the preceding H's. See Jl]| for explicit expressions for 
the G a . The important consequence of the above structure is that after eliminating the K s 
from the holonomy held equations ([81]) we obtain a very attractive set of three differential 
equations for H a : 


d 2 H n 


.dH n 


dK 2 +X dR + YHa ^dR 2 ' " dR ’ 


d 2 G n ()G r 
+ v- 


(90) 


where A", Y, /i and v are functions only of Ax and Ax. 

These three equations are linear with the only difference among them being in the homo¬ 
geneous terms: the first of them, the equation for Hi, is homogeneous with the coefficients 
depending only on Ai (a function of R), while for the next pair the inhomogeneous terms 
are driven by the solutions of the previous ones. From this structure, if the homogeneous 
equation could be solved, the remaining ones could be solved by quadratures. 

To summarize, we have so far used the information contained in two (of the three) 
Bianchi identities, one of the (A, ih)-equation, and all the held equations. The third Bianchi 
identity has not yet played any role. Assuming that we could solve (90) for H\ as a function 
(functional) of Ai we would have all the H a = H a ( Ai,Ai) and their complex conjugates. 
These expressions, which would be explicit functions of and A 1; could then be substituted 
in the third Bianchi identity. The resulting equation would involve only A, A, and the data. 
The single resulting equation, involving only Ax and A 1; i.e., 


d(H 3 /5 + G 3 ) - m 3 + 2 (H 2 /8 + G 2 )(H 3 - 1) - 2 H 2 (H 3 /5 + G 3 ) 

= -Su s + 2a B H 2 , 


(91) 


would be our sought for equation for the determination of the characteristic surfaces and 
therefore the conformal metric of the Einstein spacetime. 

The equation, however, can be manipulated, from the leading terms (but not explicitly) 
of the remaining A (H) -relations, into a simpler and more attractive equation that is referred 
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to as the light cone cut equation (LCCE) whose solution directly yields the light cone cut 
function. 

Using the solution to equation © in the form 

H\ = ^Ai + 0(A 2 ), (92) 

we have from substitution into equation (pl|) 

5 3 Ai = —45a s + higher order terms, (93) 


where a B is free data given on X + . 
A(i/)-relations we obtain, 


Manipulating this and using the equation fl85|) of the 


5 2 A 0 = Q 2 a B + 0 2 (T s + higher order terms, 
or, from the definition d 2 Z t o = Ao, 

5 2 d 2 Z = d 2 a B + 5 2 cr s + higher order terms 


(94) 


(95) 


as our final equation for the light cone cut function Z. Although the higher order terms have 
not been worked out explicitly, they can be obtained via a perturbation scheme (see next 
Section). The solutions to equation (|h|) yield the conformal structure of asymptotically flat 
vacuum spacetime. 

In the above analysis it was assumed that the solution to the equation for H\ (A), namely, 


d 2 H, 

dR 2 


A-f| + rai = 0 , 


(96) 


was known. Since it is a linear second order o.d.e. solutions for arbitrary A (R) must exist, 
though it is not clear-nor is it likcly-that we will be able to solve it explicitly. It will however, 
always be of the form 

Hi = ^Ai + 0(A 2 ) (97) 


with the higher order terms computable perturbatively. 
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Since the LCCE, (|95| ) , determines the conformal structure of vacuum spacetimes, the 


only remaining quantity to be determined is the conformal factor needed to convert the 
conformal metric into a vacuum metric. In the 9 l coordinate system the relation between 
the conformal metric g and the full metric g is given by 

g ij = tfgv (98) 

It can be shown |l[ that the conformal factor can be expressed as a simple function of H i 
and II \. namely, 

VI 2 = (l + h 1 )(l + h 1 ) +Hjh. (99) 


[As an aside we remark that from equation (|96|) , it is possible to show that the conformal 


factor (|99|) satisfies the equation 


d 2 n 

ZR2 


= Q(A)fi, 


where Q = Q( A) is a simple known function of A^ The above equation is known in the 
literature [|13|] as the Einstein bundle equation.] 


IV. AN ITERATIVE SCHEME. 

Much of the understanding of the structure and meaning of our final equations has come 
from the application of an iterative scheme to the formal theory and looking at the leading 
behavior. In this section, we present the results obtained in linear order. 

We expand all quantities in powers of a small parameter e which measures the deviation 
from flatness. From the assumption that e enters as a multiplicative factor of the Bondi 
shear, i.e. via cog, it becomes clear that the expansions have the form 


( 0 ) ( 1 ) . ( 2 ) 

Z — Z + e Z + e Z 4-..., 

i mi mi 2 ( 2 )* 

e a = e a + ee a + e e a + '"> 


U) 2 (2) 


( 3 ) 


A — eA Te A + e A..., 


( 100 ) 

( 101 ) 

( 102 ) 
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(103) 

(104) 


H 

h 


(i) 

e H 


( 2 ) 


( 3 ) 


e H +e 6 H ..., 

( 3 ) 


„ ( 2 ) „ . 

£ h + e h ~\~ _ 


There is no zeroth order contribution to A or the holonomy operator since these quantities 
are zero for Minkowski space. Furthermore, a direct calculation (using equations (|86|), (|88D, 
and (|92|)) shows that h begins at second order. 

Since A starts with order one in the perturbation expansion, it follows from the relation¬ 
ship between Z and A that ^ satisfies 

( 0 ) 

d 2 Z a = 0, (105) 

whose solution is the Minkowski space light cone cut function 


with 


( 0 ) 


(C, C) = ^(1 + CC, C + C,*(C - 0, -1 + CO; p=(i + CO- 


A. Linearized Gravity. 

We begin the linearization with the three holonomy Bianchi identities, which to first 
order are: 

(i) 1 _ (i) 

H 2 = -QH 1 , (106) 

(i) _ (i) 

H 3 = -m 2 , (107) 

m 3 = dw B , (los) 

where the K s have disappeared since they have no first order contribution. The next step in 

the process of elimination would naturally be to solve the holonomy field equations obtaining 
the H 1 s in terms of the A. However, to first (and second order) we can circumvent this step 
and consider the first of the (A, i/)-equations directly. It leads to an expression for H \, 
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which with equations ( |106| ) and (|107|) , yields the hrst order H a that automatically satisfy 
the held equations ( |90| ) to hrst order. (For third and higher order we would have to solve 
the held equations before using the A(iJ)-relations for the simplification procedure.) 


We thus have 


(i) 1 (i) 

Hl= 2 Al ’ 


(1) 1-(1) 

H 2 =-J3A U 


H 3 = -h 2 a\. 


(109) 

( 110 ) 
(111) 


Our hnal equation for the cut function is the linearized third Bianchi identity, which, after 
substituting for i/ 3 , is 


5% = —4<3(f s . 


( 112 ) 


(Although we did not have to solve the holonomy held equations, it is a straightforward 
but slightly tedious process to verify that the expressions (|109|)-(|TTT|) for the H’s do satisfy 


the held equations. A minor subtlety in the calculation is that the linearization of the held 
equations yields identically vanishing expressions, so that one must look at the equations 
at their hrst non-vanishing order. Performing the calculation, nevertheless, serves as a 
consistency check between the different sets of equations.) 

At this point we have used all the equations except the three remaining A(i/)-relations, 


A- = isAi, 

A + = -l8A„ 

Aq — A s + - A, — —DSAj. 


(113) 

(114) 

(115) 


These relations can be used to write (|112|) in the more symmetric form as the linearized 
version of the light cone cut equation: 


9 2 3 2 Z = d 2 a B (Z) + d 2 a B (Z), 


(116) 
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which is the light cone cut equation accurate to first order. In other words, this equation 
is equivalent to the the linearized conformal Einstein equations. (The above equation for 
linear theory has been derived from the vanishing of the Bach tensor, by Lionel Mason [|T4|j.) 

In linear theory (see(99))the conformal factor equals 1, and therefore in this approxima¬ 
tion we have not only the conformal Einstein equations but the Einstein equations them¬ 
selves. 

Equation ( |116|) , is to be thought of as an equation for Z, whose solution can be written 


as a sphere integral using the Green function G of the operator 5 2 9 2 . See Appendix C of 
for a derivation of the Green function of this operator. In other words we can write 


( 0 ) 


( 0 ) 


Z(x, c,<) = / G(C,C;tJ,t))[3 a,(Z + B a B (Z (x, 17 , rf))]dS v , 


(117) 


with dS v the sphere volume element in (r),rj) coordinates, as the general solution to the 
asymptotically flat linearized vacuum Einstein equations. In principle, though we have not 
yet done so in practice, the higher order terms could be calculated successively by similar 
integrals but with the integrands depending on terms to lower order in the perturbation 
scheme. 


V. SUMMARY AND CONCLUSIONS. 

We have obtained three sets of coupled equations, the holonomy Bianchi Identities, 
the A(i/)-relations and the “held equations” for the holonomy operator H and the light 
cone cut function Z, which are equivalent to the full vacuum Einstein equations. These 
equations, which already have built into them the free choice of Bondi data a B (u,(,C), can 
be manipulated and simplified (in structure) to one (complicated) equation, the LCCE, and 
a simple one for the conformal factor. On analysis, they yield perturbatively, a D’Adhemar- 
like formulation of general relativity. 

On the negative side these equations are quite unusual and are based on unfamiliar 
ideas and variables, and unfortunately are quite complicated. They nevertheless have some 
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features of considerable attractiveness: new insights often can be gained from the use of 
new variables; the perturbative solutions, from given data, are essentially unique; as the 
perturbation calculation proceeds to higher order the formalism yields the corrected light 
cone structure from the preceding order, in contrast to the usual perturbation theory which 
uses, at all orders, the Minkowski light cone structure. 

As a final comment, we remark that the work reported on here is the direct antecedent 
of another approach to general relativity (the Null Surface Theory of GR, to be reported 
elsewhere) that is based solely on families of characteristic surfaces as the basic variable of 
the theory-without any mention of fields or holonomies. The new view, which has certain 
similarities to the present work and will probably yield considerable simplifications over it, 
could not have been developed without the current approach. 

Acknowledgements: This work was written under a collaboration grant from the 
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